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Abstract] Sarcosaphagous flies are frequently found on dead bodies shortly after death. Species within this family differ in their 


developmental times[] so an accurate identification of every species is necessary for the correct estimation of the post mortem 


interval] PMI[]. Identification of specimens is traditionally performed using their morphological features. The morphological 


similarity of sarcosaphagous flies especially their eggs[] larva and pupae[] poses a challenge for forensic entomologists. Therefore 


a molecular method was established for species identification. In this study[] A 278 base pair region of mitochondrial DNA 











O Fabricius] Lucilia sericatd] Meigen[] 








O mtDNA[] coding for cytochrome oxidase subunit ond] CO [ [] was investigated for identification of the following forensically 
important species of sarcosaphagous flies from western Chinal] including Ophyra capensi{] Wiedemann[]] Chrysomya megacephala 





Lucilia cuprind] Wiedemann[[] and Boettcherisca peregrind] Robineau-Desvoidyl]. The 


results indicated that the CO | region sequenced allowed identification of major species[] providing separation of congeneric 








species with high support. In contrast[] 





and L. cuprina groups because of low sequence divergence. 


the data could not distinguish taxa from the same species group[] i.e. the L. sericata 


Key words|] Sarcosaphagous flies[] species identification[] mitochondrial DNA[] cytochrome oxidase | 


At the end of the 19th century[] the French 
entomologist Megnin and Benecke used sarcosaphagous 
flies that had developed on dead horses lying in a field 
to gather forensic informatio] Amendt et al .[] 19980. 
Todayl] forensic entomology[] including the study of 
ecology and development of sarcosaphagous flies 
colonizing a corpse[] is commonly used in investigations 
of death] Leclereq{] 1978]. In practice[] the accuracy 
of this method is about 1 day for a 1-month-old corpse 
and 1 week for a 6-month-old corps] Gaudry et al .[ 
20011. 

The first step[] however[] is to identify accurately 
sarcosaphagous flies[] which is conventionally done by 
using flies’ morphological features. The morphological 
similarity poses a great challenge for forensic 
entomologists[] especially in identifying larvae. The 
minuscule morphological differences between species of 
immature larvae are sometimes extremely hard to 
determine and requires an expert with specialized 
taxonomic knowledge in the field to assure correct 
identification [] Catts and Goff] 1992[]. Although 
identification keys are available[] only a few experts are 
able to identify the larvae of sarcosaphagous flies to 
species level. Up to now[] all eggs and pupa and the 





first to third instar larvae of sarcosaphagous flies could 
not be identified by morphologyical method[] Wallman 
and Donnellan{] 2001[]. Under these circumstances[] 
species identification based on genetic examination is an 
option. 

Such morphological challenges make a molecular 
approach to the identification of these species desirable. 
Recent molecular studies have demonstrated the 
applicability of mtDNA sequencing to the identification 
of important Old and New World species of carrion- 
breeding blowflies[] Catts and Goff{] 1992]. When 
using mtDNA sequences for species identification[] one 
must consider the possibility of a single sequence may 
not be representative of the entire species. However[] 
previous studies of mtDNA restriction site variation 
within a wide variety of insect species indicate that the 
extent of variation within species is almost always small 
compared to differences between species[] Sperling[] 
1993[]. For this study of divergence between species|] 
people deliberately chose sequences that are known to 
show relatively low levels of variation within populations 

O Simon[] 1991[]. The recent study indicated that 
analysis of the partial sequences focused on sections of 
the cytochrome oxidase b subunit one and two] CO [| 
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and COII [] encoding region of mtDNA revealed 
abundant phylogenetically informative nucleotide 
substitutions that could identify sarcosaphagous flies at 
species to species group level. Our studies involved 
sequencing to examine the potential of DNA for use in 
the identification of blowfly in forensic investigation. 
The sequence variation within the CO | gene of the 
mitochondrial genome of several species was selected for 
examination. Analyses were based on a 278 base pair 
region of the gene coding for cytochrome oxidase subunit 


ond] CO I D. 
1 MATERIALS AND METHODS 


1.1 Insects 
Four adults of Chrysomya megacephala 
[ Fabricius] 2 adults of Lucilia sericata] Meigen[] and 
their eggs[] 2 adults of Boettcherisca peregrina 
O Robineau Desvoidy[] and its 1 larva and 1 pupal] 4 
adults of Lucilia cuprinal] Wiedemann{[] 2 adults of 
Ophyra capensis || Widemann[] were obtained from 
Huhhot district SAT] 39°35’ S[] 110°46’ ELL] western 
China. In addition[] specimens from Chengdu and 
Dunhuang district of western China were included. 
These were Chrysomya megacephala [] 4 adults [O 
Lucilia sericata] 2 adults[[] Boettcherisca peregrinal] 2 
adults[[] Musca domestica vicina[] 2 adults[] were 
obtained from Chengdu district] SA] 30°05’S{] 102°54’ 
EQ) and Boettcherisca peregrina |] 1 adult[] from 
Dunhuang district] SA] 37°56’ SQ 92°07'EL]. 
1.2 DNA extraction 
DNA was extracted using Chelex method] Walsh 
and Petzger[] 1991[]. DNA was quantified through 
Primate-specific alpha-satellite probe assay[] Waye et 
al .[ 19890. 
1.3 PCR primers 
A portion of the CO | gene was amplified and 
sequenced using primers C1-J-2499] 5’-CAG CTA CTT 
TAT GAG CTT TAG G-3'[[] sense[] and C1-N-2800 
O 5'-CAT TTC AAG CTG TGT AAG CAT C-3’[[] anti- 
sensel]] Sperling et al .[] 1994[] Wells and Sperling{] 
1999[]. 
1.4 PCR conditions 
The PCR reaction volume was 37.5 „UO 
containing 20 — 40 ng DNA[] 6 pL dNTP 1 mmol/ 
mL[]]1.5 U Taq polymerase[] 3.75 uL 10 x buffer 
O Mg** 1.5 mmol/L 0.3 pL primer{] 50 nmol/mL{. 
PCR amplifications were performed in a 
thermocycler[] Perkin-Elmer 9600[[] with initiative 
denaturing for 1 min at 94°C] followed by 38 cycles of 
94°C for 1 min[] 48°C for 1 min and 72°C for 2 min. 
1.5 Sequncing 
After purification of the PCR products with QIA- 
quick columns cycle sequencing was performed on both 
forward and reverse strands using ABI PRISM Big Dye 


Terminator Cycle Sequencing Ready Reaction Kit 
O Applied Biosystems []. Removal of excess dye- 
deoxyterminators primers and buffer was accomplished 
by DYE-EX spincolumnd] Qiagen[]. 
1.6 Phylogenetic analysis 
Phylogenetic analysis is the study of the 
relationships between different groups of organisms[] and 
is thus based on similarity and difference in characters 
choser{] Hillis and Moritz[] 1990[]. As the sequences 
were protein coding and did not contain any insertions 
or deletions[] alignments were carried out using the 
multiple-alignment program DNAMAN[] version 4.0 
analysesL] Lynnon BioSoft{] Copyright 1994-98[]. We 
tested whether the sequences were of mitochondrial 
origin or represented paralogous sequences resident in 
the nucleus{] Hillis et al .[] 1996[] in two ways using 
MEGA[] Kumar et al.[] 1993[]. Sequences were 
translated with mitochondrial genetic code and[] since 
nuclear paralogues can lose their coding function[] the 
predicted amino acid sequences were inspected for 
inappropriate stop codons. We also compared the base 
composition of the individual sequences since nuclear 
paralogues can have divergent base compositions relative 
to mitochondrial genes. A neighbour-joining tree using 
the Tamura and Nei] Tamura and Meil] 1993[] model of 
nucleotide substitution was constructed using the MEGA 
packag¢] Kumar et al .[] 1993[]. Maximum parsimony 
analysis was also conducted to calculate the tree based 
on the least number of evolutionary steps. Five hundred 
bootstrap replications were performed. Bootstrap values 
are calculated by repeated random sampling of the data 
to provide an indication of the confidence limited for a 
particular group. 


2 RESULTS 


2.1 Morphological identification of specimens 
Eggs obtained from abdomen of adults of Lucilia 
sericata were definitely known for their species identity . 
But one pupa and larva] first instar{] from a rabbit body 
[ in the experiment [] couldn’t be identified by 
morphological method] Hu and Min[] 20000. All adults 
of specimens were identified easily by morphological 
method except the two most closely related species[] i. 
e. Lucilia sericata and Lucilia cuprina . 
2.2 COT sequence divergence 
Six species were sequenced. The sequences have 
been deposited in GenBank with accession number 
AY818091 — 818125. 
2.3 COT sequence analysis 
The neighbor-joining] NJ[ tree based on the CO 
I sequence data using Kimura’s two-parameter model 
was constructed|] Fig. 1[]. Invariant positions in the 
sequence were removed and the remaining variant 
positions numbered relative to the original sequence to 
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indicate the number of substitutions supporting the 
divergence of the taxa. The neighbour-joining tree was 
identical to the tree obtained using maximum parsimony 
method. Three distinct congeneric clusters were formed 
based on the sequence data. High bootstrap values 
support the three nodes. Bootstrap values provide an 
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0.01 


Fig. 1 


indication of the percentage support for a grouping by 
randomly resampling the data. 
2.4 Phylogenetic analysis 


A total of 278 aligned sites for the 29 CO Į 
sequences were included in the analysed] Fig. 1{]. 
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The neighbor-joining] NJL] tree based on the CO | sequence data using Kimura’s two-parameter model 


The bar indicates 0.01 substitutions per site. 


All specimens were grouped into 4 clusters with 
high bootstrap support. The first and second clusters 
were formed by 2 groups respectively[] and the third and 
fourth by 1 group respectively. At the clusters level{] 
there was 98% support between the first and the second 
clusters. At the group level within the first cluster[] 
there was 28% support between the first and second 
groups. At the species level[] the first group of the first 
cluster was formed by 3 subgroups of specimens of 
Lucilia sericata and Lucilia cuprina. Within the first 
subgroupl] the first 3 specimens of Lucilia sericata 
shared 99% support[] whereas the fourth specimen 
deviated a little as showed as 97% . Within the second 
subgroup containing 2 species[] i.e. Lucilia sericata 
and Lucilia cuprina[] the support between species is 


surprisingly greater than the support within species as 
shown as 71% vs 99%. Within the third subgroup 
there was 99% bootstrap support between the two 
specimens of Lucilia cuprina from Huhhot. The 
bootstrap support among all three sub-groups was 75% . 
In the second group of the first cluster[] there were 8 
of Chrysomya megacephala with 
approximately 52% support among them. 

In the second cluster{] there were 2 groups[] with 
99% support between them. The support within the 
first group was as high as 99% whereas the support 
within the second group of adult specimens was 99% [] 
and the support between adult specimens as a whole and 
that of larva was 82% . 

In the third and fourth clusters[] the support was 


specimens 


30 CAI Ji-Feng et al .[] The availability of mitochondrial DNA cytochrome oxidase | gene for the distinction 





of forensically important flies in China 383 
99% respectively at the species level. 2.8%[] 2.86 %[] 3.59%[] 0.37%[] and O 
Table 1 shows that the level of CO | gene respectively. While the maximum and minimum level of 
nucleotide divergence between species groups within divergence between the Chrysomya megacephala and 
sarcosaphagous flies taxa. A total of 278 aligned sites Lucilia sericata was 12.6% and 6.4% respectively 
for the 29 CO | sequences were included in the  [] average 8.9%[[] that of Lucilia sericata and Lucilia 
phylogenetic analyses by Kimura 2-parameter distance cuprina was 8.5% and 0.4% respectively[] average 


method. Uncorrected percentage sequence divergence 5.1%] that of Chrysomya megacephala and 
among all taxa ranged from 0 to 17.8%. The overall Boettcherisca peregrina was 13.4% and 11.3% 
average was 9.9% [] and the average level of divergence respectively[] average 12.7%[]] that of Chrysomya 
within the Chrysomya megacephalal] Lucilia sericata[] megacephala and Musca domestica vicina was 17.0% 
Lucilia cuprina[] Boettcherisca peregrina] Ophyra and 15.2% respectively] average 15.4%[. 

capensis[] Musca domestica vicina group was 0.29% [] 


Table 1 Percent uncorrected sequence divergence between sarcosaphagous flies taxa for the CO] gene’ 


1 2 3 4 5 6 7 8 9 10 Il 12 13 14 15 16 17 18 19 20 21 2 23 24 25 25 27 28 29 







































































Oi - 

O20 1.50 - 

O30 1.50 0.00 - 

O 4 1.50 0.00 0.00 - 

O50 1.50 0.00 0.00 0.00 - 

O 6) 1.50 0.00 0.00 0.00 0.00 - 

O 70 1.50 0.00 0.00 0.00 0.00 0.00 - 

O 80 1.50 0.00 0.00 0.00 0.00 0.00 0.00 - 

0 90 9.20 8.00 8.00 8.00 8.00 8.00 8.00 8.00 - 

0 100 7.60 6.40 6.40 6.40 6.40 6.40 6.40 6.40 3.00 - 

O 110) 12.6 11.3 11.3 11.3 11.3 11.3 11.3 11.3 8.90 6.10 - 

O 120 9.20 8.00 8.00 8.00 8.00 8.00 8.00 8.00 0.00 3.00 8.90 - 

0 130 12.6 11.3 11.3 11.3 11.3 11.3 11.3 11.3 8.90 6.10 0.00 8.90 - 

O 140 9.20 8.00 8.00 8.00 8.00 8.00 8.00 8.00 0.00 3.00 8.90 0.00 8.90 - 

[0 150 8.00 6.80 6.80 6.80 6.80 6.80 6.80 6.80 5.30 3.70 7.30 5.30 7.30 5.30 - 

0 160 12.2 10.9 10.9 10.9 10.9 10.9 10.9 10.9 8.50 5.70 0.40 8.50 0.40 8.50 6.80 - 

O 170 7.60 7.02 7.02 7.20 7.20 7.20 7.20 7.20 5.70 4.10 7.70 5.70 7.70 5.70 0.40 7.30 - 

0 180 12.2 10.9 10.9 10.9 10.9 10.9 10.9 10.9 8.50 5.70 0.40 8.50 0.40 8.50 6.80 0.00 7.30 - 

0 190 13.0 11.3 11.3 11.3 11.3 11.3 11.3 11.3 11.7 11.3 10.9 11.7 10.9 11.7 11.7 10.9 12.1 10.9 - 

0 200 13.4 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 11.7 13.1 11.3 13.1 3 13.0 13.1 13.5 13.1 4.50 - 

0 210 12.6 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 11.3 11.7 11.3 11.7 3 11.3 13.0 11.7 11.3 8.40 6.80 - 

0 220 13.4 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 11.7 13.1 11.3 13.1 11.3 13.0 13.1 13.5 13.1 4.50 0.00 6.80 - 

O 230 13.4 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 11.7 13.1 11.3 13.1 3 13.0 13.1 13.5 13.1 4.50 0.00 6.80 0.00 - 

O 240 12.6 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 11.3 11.7 11.3 11.7 3 11.3 11.3 11.7 11.3 8.40 6.80 0.00 6.80 6.80 - 

0 250 12.6 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 11.3 11.7 11.3 11.7 11.3 11.3 11.3 11.7 11.3 8.40 6.80 0.00 6.80 6.80 0.00 - 

0 260 12.6 11.7 11.7 11.7 11.7 11.7 11.7 11.7 11.3 10.4 10.8 11.3 10.8 11.3 8.00 10.8 8.40 10.8 12.1 14.3 15.1 14.3 14.3 15.1 15.1 - 

0 270 13.1 12.2 12.2 12.2 12.2 12.2 12.2 12.2 11.7 11.7 12.1 11.7 12.1 11.7 9.20 12.1 9.60 12.1 12.5 13.8 14.7 13.8 13.8 14.7 14.7 11.0 - 

0 280 17.0 15.2 15.2 15.2 15.2 15.2 15.2 15.2 16.0 16.0 16.9 16.0 16.9 16.0 14.2 16.4 14.7 16.4 16.0 17.4 17.8 17.4 17.4 17.8 17.8 12.5 12.1 - 

0 290 17.0 15.2 15.2 15.2 15.2 15.2 15.2 15.2 16.0 16.0 16.9 16.0 16.9 16.0 14.2 16.4 14.7 16.4 16.0 17.4 17.8 17.4 17.4 17.8 17.8 12.5 12.1 0.00 - 
For species that diverged from all others by a relatively large percentage and where within-species variation was small[] the minimum sequence divergence percentages are given that 
distinguish any individual within that species from any other species. Species/specimen codes{] 1 - 4[] Chrysomya megacephald] Huhhot{[] 5 - 81] Chrysomya megacephald] Chengdu 
90100 Lucilia sericatd] Chengdu[[] 110120 Lucilia sericatd] Huhbot{]] 130140 Lucilia sericata eg] Hubhot{]] 15 - 180 Lucilia cuprind] Hubhot[[] 19{] Boeticherisca peregrina larva 








O Huhhof l] 200 Boettcherisca peregrina pupa] Huhhot{]] 210) Boettcherisca peregrinal] Dunhuang{[] 220230 Boettcherisca peregrinal] Hubhot{]] 240250 Boettcherisca peregrina larva 
O Chengdu 26270 Ophyra capensid] Huhhoi] 28290 Musca domestica vicind] Chengdu. 


the difficulty in differentiating the closely related 


3 DISCUSSION species shows Lucilia sericata and Lucilia cuprina . 

However[] the morphological method requires expertise 

In this study[] the DNA extraction and sequencing in specialized taxonomy to assure correct identification 

of all specimens of six species of sarcosaphagous flies [] Leclercq] 1978[[] whereas the technology of mtDNA 

were finished. Both specimens from Huhhot that were sequence method turns to be mature and the concerned 

stored in 75% alcohol for 3 - 6 months and from facilities are easily accessible even in ordinary 
Chengdu that were stored by air-dried condition for over laboratories . 

2 years can be extracted and sequenced. According to the neighbor-joining] NJ[ tree based 


The results indicate that the morphological method on the CO | sequence data by using Kimura’s_ two- 
is as effective as mtDNA sequence method in species parameter model [] Fig. 1 [O Calliphoridae 
identification of adult sarcosaphagous flies and can solve Sarcophagidae and Muscidae of Diptera clustered in 
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divergent branches. The high support for the congeneric 
grouping of species illustrates the potential of the CO I 
in interspecific distinction. The ability to clearly 
distinguish between these forensically prominent genera 
based on such a small region provides a strong 
indication of economy in choosing the region of the CO 
I. The difficulty we encountered in identifying the 
specimens of Lucilia sericata and Lucilia cuprina by 
mtDNA sequence method is corresponding to that 
reported by Wallman and Donnellar{] Wallman et al. [0 
20010] who encountered the difficulty in distinction 
between C. augur and C. dubia based on CO [ 
sequence data in their study[] and also the report of 
Stevens and Wall[] Stevens and Wall[] 1996[[] who 
examined the mitochondrial encoded 12S rRNA[] CO I 
and CO [| sequences of the two blowflies Lucilia 
cuprina and L. sericata originated from Hawaii. The 
two endemic Chinese species in our study were closely 
related to each other and further sequencing is desired 
to separate them. The intraspecific variation within 
Lucilia sericata and Lucilia cuprina should deserve 
further consideration . 

The results of our study about CO [ phylogenetic 
systematic data are in agreement with Wells and 
Sperling’{] Wells and Sperling[] 1999[] Wells et al .{] 
2001 [[] who demonstrated] by examining CO | 
sequences|] that the blowflies Chrysomya rufifacies and 
Chrysomya albiceps exhibited less than 1% intraspecific 
difference and about 3% interspecific one. However[] 
the results of our study are in disagreement with their 
results concerning Lucilia sericata[] Lucilia cuprina[] 
Boettcherisca peregrina groups. The difference may be 
due to geographical separation of populations[] which 
needs further study . 

Though some studies reported the differentiation of 
samples from different fields[] our study showed no 
differentiation between specimens from Dunhuang and 
Huhhot and lent little support to the differentiation 
between specimens from Chengdu and Huhhot in 
western China. So in the future[] we should sample 
more specimens from more species and from different 
geographical regions throughout the world. 

This report shows that[] with future development] 
mitochondrial nucleotide sequences can provide a means 
of identifying forensically important flies. Sequencing 
this region of mtDNA has provided data that appear to 
indicate the potential of mtDNA in development of an 
identification protocol for use in forensic entomological 
studies. Future works include the examination of more 
variable mitochondrial genes[] if necessary[] and the 
analysis of additional species of flies from other parts of 


China. 
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